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Abstract 

We investigate the potential energy curves of rare-gas dimers with various ranges and strengths 
of interparticle interactions (nuclear-electron, electron-electron, and nuclear-nuclear interactions). 
Our investigation is based on the highly accurate coupled-cluster theory associated with those inter¬ 
particle interactions. For comparison, the performance of the corresponding Hartree-Fock theory, 
second-order Mpller-Plesset perturbation theory, and density functional theory is also investigated. 
Our results reveal that when the interparticle interactions retain the long-range Coulomb tails, the 
nature of van der Waals interactions in the rare-gas dimers remains similar. By contrast, when the 
interparticle interactions are sufficiently short-range, the conventional van der Waals interactions 
in the rare-gas dimers completely disappear, yielding purely repulsive potential energy curves. 
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I. INTRODUCTION 


Van der Waals (vdW) interactions [1-5] are omnipresent in materials and biological sys¬ 
tems. These interactions are of fundamental importance in numerous fields, involving molec¬ 
ular and condensed matter physics, supramolecular chemistry, structural biology, surface 
science, and nanoscience. While vdW interactions are individually weak (e.g., compared to 
covalent bonds or electrostatic interactions between permanent charges, dipoles, etc.), they 
are collectively important in the determination of the structure, stability, and function of a 
vast variety of systems, such as the interaction between graphene layers, the self-assembly 
of functional nanomaterials, the structure of biomacromolecules (e.g., DNA, RNA, and pro¬ 
teins), and the molecular recognition of proteins 6]. 

In particular, the potential energy curve of a rare-gas dimer is predominantly determined 
by the interplay between the exchange-repulsion energy at short internuclear distances and 
the attractive vdW interaction at large internuclear distances, exhibiting a potential mini¬ 
mum (the vdW minimum) at an intermediate internuclear distance. The exchange-repulsion 
energy arises from the overlap of the electron densities of the two atoms. On the other hand, 
the vdW interaction, also known as London dispersion interaction or induced dipole-induced 
dipole interaction, arises from the Coulomb correlation of electron density fluctuations in 


the two well-separated atoms. The po 
by the Lennard-Jones (LJ) potential 


;ential energy curve can be conveniently approximated 

1 


Vlj (R) = 4e 


12 


( 1 ) 


where R is the internuclear distance, cr is the distance at which the potential is zero, and —e 
is the minimum of the potential, which is reached at R — 2 1 / 6 cr. Here the R~ 12 term models 
the exchange-repulsion energy, dominant at short internuclear distances, while the R ~ 6 term 
models the attractive vdW interaction, dominant at large internuclear distances. Whereas 
the attractive term is physically based, the repulsive term has no theoretical justification 
(i.e., chosen for computational efficiency). Note that the exchange-repulsion energy should 
decay almost exponentially with the internuclear distance. Nevertheless, due to its com¬ 
putational simplicity, the LJ potential is widely used in computer simulations even though 
more accurate potentials exist. 

However, the R~ 6 dependence of the vdW interaction may not be applied to macroscopic 
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systems like colloids and biological membranes. In these systems, the vdW interaction 
between two objects immersed in a medium is strongly influenced by the dielectric properties 
of the objects and the medium. Accordingly, the resulting vdW interaction can be very 


under certain conditions 


0 , 00 . 


different from the conventional R 6 expression 2 ), 7, [8j, and can be completely repulsive 


Severa 


ascinating phenomena have been discovered in 


12 ], 


these non-fU 6 macroscopic vdW systems 

Is it possible to create non-f?~ 6 vdW interactions between rare-gas atoms in vacuum? 
Conceptually, the types of interparticle interactions (nuclear-electron, electron-electron, and 
nuclear-nuclear interactions), traditionally given by the Coulomb interactions, should play 
a fundamental role in determining the properties of atoms and molecules. Hence, we expect 
that non-i? -6 vdW interactions can appear by tuning the effective interparticle interactions 
of rare-gas atoms in vacuum. As a proof of concept, in this work, we address how the nature 
of vdW interactions in rare-gas dimers (i.e., the simplest vdW systems) changes with varying 
interparticle interactions, using the highly accurate coupled-cluster theory associated with 
those interparticle interactions. The rest of this paper is organized as follows. In Section II, 
we describe our model systems and computational details. We compare the results obtained 
from the coupled-cluster theory with those obtained from different computational methods, 
and give our comments on the connection between this study and a popular scheme in 
density functional theory in Section III. Our conclusions are given in Section IV. 


II. MODEL SYSTEMS AND COMPUTATIONAL DETAILS 


For a system consisting of M nuclei and N electrons in the Born-Oppenheimer approx¬ 
imation (as the nuclei are much heavier than the electrons), the electronic Hamiltonian 
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H e = - 
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is the sum of the kinetic energy of electrons, the nuclear-electron attraction energy, and the 
electron-electron repulsion energy, respectively. Here Za is the atomic number of nucleus A, 
m e is the mass of an electron, — e is the charge of an electron, = |r^ — R 4 I is the distance 
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between electron i and nucleus A, r ^ = |r* — i\,-| is the distance between electrons i and 
j , and /(r) is the interparticle interaction operator with r being the interparticle distance. 
The electronic Schrodinger equation 


H e ^ e = £ e T e 


( 3 ) 


is solved for the electronic energy E e and electronic wavefunction T e , which describes the 
motion of the electrons for fixed nuclear positions. The total energy 


M M 


Efotal ~ E e + 


47re f 
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( 4 ) 


A= 1 B>A 

is obtained by adding the nuclear-nuclear repulsion energy to the electronic energy, where 
R ab = |Ra — Rb| is the distance between nuclei A and B. One can obtain E to tai as a function 
of the nuclear positions, commonly known as the potential energy curve (or surface). 

Traditionally, /(r) is given by the Coulomb interaction 1/r. However, in this work, we 


ojr)/r, which are generated by splitting the 
14], The former (the erf interaction) retains 


consider two types of /(r): exi(ujr)/r and erfc 
Coulomb interaction into two components 
the long-range Coulomb tail without the singularity at r = 0, while the latter (the erfc 
interaction) is a short-range interaction with a singularity at r = 0. Physically, 1/uo spec¬ 
ifies the distance beyond which erf(cur)/r approaches 1/r and the distance beyond which 
erfc {ujr)/r becomes insignificant (see Figure 111). Similar to the Coulomb case 15, ll6], the 
nuclear-attraction and two-electron repulsion integrals modified for the erf and erfc inter¬ 
action operators can be evaluated analytically over Gaussian basis functions 16-18], facil¬ 
itating an efficient evaluation of the integrals needed for solving Eq. (|3|) and the equations 
associated with related approximate methods (see below). In principle, other types of fir) 


can also be adopted 


HE 


22 ] 


Similar to the Coulomb case, solving Eq. (|3|) for a given /(r) is, however, extremely 
difficult even for the ground-state energy and wavefunction of a very small system, due to 
the prohibitively expensive computational cost. Practically, one searches_ for approximate 
solutions to Eq. (|3j), obtained from ab initio wavefunction methods jp, is], 1^, 23], such as the 
Hartree-Fock (HF) theory, second-order Mpllcr-Plesset perturbation theory (MP2), coupled- 
cluster theory with iterative singles and doubles (CCSD), and CCSD with a perturbative 
treatment of triple substitutions (CCSD(T)). Among them, the CCSD(T) method with a 
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sufficiently large basis set is generally expected to provide highly accurate results for a variety 
of small- to medium-sized systems. 

Alternatively, Kohn-Sham density functional theory (KS-DFT) 24], a popular method 
for the study of the ground-state properties of large systems, can also be employed. Similar 
to the Coulomb case, density functional approximations (DFAs), such as the local density 
approximation (LDA) and generalized-gradient approximations (GGAs), to the exchange- 
correlation (XC) energy functional for a given f(r) are needed in the corresponding KS-DFT 


25 


26]. Here, the LDA exchange energy functional for the erf interaction is obtained by 


subtracting the LDA exchange energy functional for the erfc interaction 27] from the LDA 


exchange energy functional for the Coulomb interaction 28], whereas the LDA correlation 


energy functional for the erfc interaction is obtained by subtracting the LDA correlation 
energy functional for the erf interaction fl front the LDA correlation energy functional for 
the Coulomb interaction |30|. In addition, as the Perdew-Burke-Ernzerhof (PBE) XC energy 
functional (i.e., a popular GGA) for the Coulomb interaction 31] and its variant for the erfc 


interaction 


32] are both available, their difference gives the PBE XC energy functional for 


the erf interaction. 

To illustrate how the nature of vdW interactions in rare-gas dimers changes with vary¬ 
ing interparticle interactions, we calculate the potential energy curves of the He-He dimer 
associated with the interparticle interactions erf (ojr)/r ( u> = oo, 10.00, 2.00, 1.70, 1.40, and 
1.10 bohr -1 ) and erfc (ujr)/r (uj = 0.00, 0.10, 0.20, 0.25, 0.30, and 0.40 bohr -1 ), using the 


corresponding CCSD(T), CCSD, MP2, HF, and KS-DF 


energy functionals for the associated interactions 
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1 employing the PBE and LDA XC 


14, 


All calculations are performed with a development^ 


are computed using a large aug-cc-pVQZ basis set |34] wit 


grid, consisting of 250 Eulcr-Maclaurin radial grid points 


points 


version of Q-Chem 4.0 33]- Results 
r a high-quality EML(250,590) 
and 590 Lebedev angular grid 


. The counterpoise correction 37] is employed to reduce basis set superposition 


error (BSSE). 


III. RESULTS AND DISCUSSION 

The potential energy curves of the He-He dimer associated with the long-range interpar¬ 
ticle interactions erf (cor)/r, calculated using the corresponding CCSD(T), are presented in 
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Figure [ 2 I. Similar to the Coulomb case (i.e., the oo = oo case of the erf interaction), all the 
potential energy curves resemble the LJ potentials. For a smaller oo, the strength of the erf 
interaction is weaker. Consequently, the electrons are more loosely bound to the nucleus, 
and the atoms are more polarizable, yielding larger values of a and e, respectively (see Eq. 
(IT]) ) 4]. Owing to the long-range nature of the erf interaction, the attractive vdW interac¬ 
tion is shown to have the [erf(u;-R)] 2 fW 6 asymptote (essentially retaining the R~ G asymptote 
of conventional vdW interactions) at sufficiently large internuclear distances R, based on a 
second-order perturbation theory (see the Appendix). 

In comparison with the highly accurate CCSD(T) results, the He-He potential energy 


ated using the corresponding CCSD, MP2, 
As shown, CCSD performs similarly to 


curves associated with the erf interactions, calcu 
HF, PBE, and LDA are presented in Figure |3 
CCSD(T), and slightly outperforms MP2. Besides, CCSD(T), CCSD, and MP2 exhibit 
the correct R~ G vdW asymptotes. By contrast, due to the lack of electron correlation, 
HF completely fails to describe the attractive vdW interactions, yielding purely repulsive 
potential energy curves for all the oo values studied. Within the framework of KS-DFT, 


PBE consistently outperforms LDA. However, in view of the large errors associated wit 


1 the 


- 6 > 


38 


vdW minima and the incorrect vdW asymptotes (decaying much faster than R 
LDA, PBE, and possibly other semilocal density functionals 40] cannot accurately describe 
long-range vdW interactions 41], wherein a fully nonlocal XC energy functional should be 


essential 


25 


26 


On the other hand, the potential energy curves of the He-He dimer associated with 
the short-range interparticle interactions erfc(o;r)/r, calculated using the corresponding 
CCSD(T), are shown in Figured. In contrast to the Coulomb case (i.e., the oo = 0 case of the 
erfc interaction), the potential energy curves show strong ^-dependence. It resembles the 
LJ potential only for a vanishingly small oo, displays a metastable state for an intermediate 
oo (around 0.25 bohr -1 or smaller), and becomes purely repulsive for a oo larger than 0.30 
bohr -1 (see the Appendix). 

For comparison, the He-He potential energy curves associated with the erfc interactions, 


y. 


calculated using the corresponding CCSD, MP2, HF, PBE, and LDA are shown in Figure 
With the increase of oo, the potential energy curves obtained from all the computational 
methods become very similar. As would be expected on physical grounds, semilocal density 


functionals can be surprisingly accurate for short-range XC effects 


40]. PBE is shown to 
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consistently perform better than LDA. Besides, due to the dominance of exchange-repulsion 
energy for a sufficiently large u>, even the HF theory can be reliably accurate. 

Similar to the Coulomb case, the overall trends of LDA and PBE are opposite to those 
of HF and MP2, implying that a combination of the HF exchange, MP2 correlation, and 
DFAs (e.g., LDA or GGAs) in KS-DFT (i.e., hybrid DFT Q or double-hybrid DFT 431) 
may achieve a more favorable balance between cost and performance than CCSD(T) for the 
vdW interactions in large rare-gas dimers under the erf and erfc interactions. 

In addition, we calculate the potential energy curves of the He-Ne and Ne-Ne dimers 
associated with the erf and erfc interactions, using the corresponding CCSD(T), as shown 


in Figures [b| to 0. For the Coulomb case, the values of a for the He-Ne and Ne-Ne dimers 
are larger than that for the He-He dimer. Nevertheless, similar trends are also found for the 
potential energy curves of the He-Ne, Ne-Ne, and possibly other rare-gas dimers. 


cu- 


44] 


To test the transferability of the above observed trends for other vdW systems, we ca’ 
late the potential energy curves for the lowest triplet states of H 2 (a simple vdW system) 
associated with the erf and erfc interactions, using the correspon ding CCSD (i.e., an exact 
theory for any two-electron system). As shown in Figures 


10 


and 


11 


, the major features of 


the potential energy curves remain very similar to those found for rare-gas dimers. 

Here we comment on the connection between this study and long-range corrected (LG) 


hybrid functionals for systems with Coulomb interactions 45 56]. These functionals model 


(i.e., orbital-c 
Figures ij| and 


the short-range interaction (e.g., the erfc interaction) by a DFA in KS-DFT and the comple¬ 
mentary long-range interaction (e.g., the erf interaction) by HF exchange or a fully nonlocal 
dependent) XC energy component from ab initio wavefunction methods. In 
5|, compared to the highly accurate CCSD(T) results, LDA and PBE perform 
reasonably well for sufficiently short-range interparticle interactions, whereas they perform 
poorly for long-range interparticle interactions. Accordingly, our findings are also in sup¬ 
port of the key feature of the LG hybrid functionals for systems with Coulomb interactions, 
which have recently been found to provide supreme performance for a very wide range of 
ications { 57 , 58], especially for problems related to the asymptote of the XC potential 


66], self-interaction errors [67], [68], fundamental gaps 69r82], and charge- 


tions 


83 


Besides, empirical atom-atom dispersion potentials 


correlation energy 


(m, 55, 


56 


ransfer excita- 


90-192] or MP2 


4311531 .193H95] can be added to the KS-DFT energy in order to improve 


the description of noncovalent interactions (e.g., vdW interactions). Alternatively, KS-DFT 
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may also be combined with symmetry-adapted perturbation theory (SAPT) 96^ 104] to 


yield accurate results for noncovalent interactions 


105- Jill]. 


In addition, to properly de¬ 


scribe strong static correlation, it could be essen 


25 


ial t o de velop a combined LC hybrid scheme 


112 


114] for small- to medium-sizec 


with random phase approximations (RPAs) 
or with thermally-assisted-occupation density functional theory (TAO-DFT) 
large-sized systems. 


115 


syst ems 
117 ] for 


IV. CONCLUSIONS 

In conclusion, we have developed a comprehensive understanding of the physics involved 
in controlling the vdW interactions in rare-gas dimers. Specifically, we have examined the 
potential energy curves of the rare-gas dimers associated with a variety of interparticle 
interactions, using the highly accurate CCSD(T) method as well as other computational 
methods. The long-range interparticle interactions are shown to be essential for retaining the 
main features of conventional vdW interactions, which cannot be properly described by LDA, 
PBE, and possibly other semilocal density functionals in KS-DFT, but can be accurately 
described by MP2, CCSD, and possibly other fully nonlocal XC energy components from ab 
initio wavefunction methods. On the other hand, the nature of vdW interactions is shown 
to change drastically with the short-range interparticle interactions, wherein LDA, PBE, 
and possibly other semilocal density functionals in KS-DFT perform reasonably well for 
sufficiently short-range interparticle interactions (e.g., erfc(cur)/r with u = 0.30 bohr -1 or 
larger). Therefore, our findings also support the main feature of the LC hybrid functionals 
for systems with Coulomb interactions. Although only the vdW interactions in rare-gas 
dimers and the triplet H 2 molecule are studied and discussed in this work, our conclusion 
may remain appropriate for other vdW-dominated systems. 
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APPENDIX: ASYMPTOTE OF THE INTERACTION ENERGY CURVE BE¬ 
TWEEN TWO WELL-SEPARATED RARE-GAS ATOMS ASSOCIATED WITH 
THE LONG-RANGE (ERF) INTERPARTICLE INTERACTIONS 


B), 


Similar to the derivation for the Coulomb case (e.g., see Chapter 3 of Ref. |5fl), we de¬ 
rive an analytical expression for the asymptote of the interaction energy curve between two 
well-separated rare-gas atoms associated with the long-range interparticle (nuclear-electron, 
electron-electron, and nuclear-nuclear) interaction operator f(r): erf(o;r)/r (the erf inter¬ 
action), based on a second-order perturbation theory [1181 ]. Since in the Coulomb case, a 
rare-gas atom has no permanent multipole moments in its nondegenerate ground state [3], 
presumably this remains correct for the erf interaction with a sufficiently large u or for the 
erfc interaction [/(r): erfc (ojr)/r\ with a sufficiently small oj. Note also that the finite speed 


of propagation of electromagnetic signals is not taken into account in our derivation 


brevity, the Einstein summation convention 


B- 


For 


119!] is adopted here. Based on this convention, 


when an index variable appears twice in a term, it implies a summation of that term over 
all possible values of the index. 

Consider a rare-gas atom A, composed of a nucleus situated at r a=0 and Na electrons 
situated at r Q ( a = 1,2, ...,Na) with respect to the nucleus of A. The electric potential at 
a point r, due to the charge distribution, is 


V A (r) = 


47re ( 


N A 


r «i), 


(5) 


a=0 


where e„ =0 = iV^e is the nuclear charge of A, and et = —e (a = 1 , 2 ,..., Na) is the charge 
of an electron. The Taylor series expansion of V^r) around the nucleus of A gives 


V A (r) = 


1 


4ne 0 

1 

47re n 


e aU<* v */(r) + 2j C>'VoCo Y ; V ; /(r 


tot 


- lit'Vi + QgViV,- + • • • 


f(r), 


( 6 ) 


where the first term is from an electric monopole ef ot = ^ e^, the second term is from an 

a 

electric dipole, whose i th Cartesian component nf = Y2 e a r ion the third term is from an 

a 

electric quadrupole source, Qfj = ^ J2et r ia r ja, an d so on. 

a 

Consider a second rare-gas atom B, composed of a nucleus situated at r i a =0 and Nb 
electrons situated at rg (/3 = 1,2, with respect to the nucleus of B. Let R be the 
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separation distance vector pointing from the nucleus of A towards the nucleus of B. The 
interaction energy between atoms A and B is 


N b 


u A B = J2 e P VA ( R+r ^’ 

0=0 


(7) 


where e^ =0 = iV^e is the nuclear charge of B, and uf = — e (/3 = 1, 2, Nb) is the charge 
of an electron. The Taylor series expansion of Va(R + r^) around the nucleus of B gives 


Va(R + r/j) = Va(R) + r i/3 ViV A (R) + —r^r^ViVjV A (R) 


1 + r i/3 Vi + ^rjoViVj + 


Substituting Eq. (JH]) and Eq. (JSJ) into Eq. (JTD produces 

1 


V a (R). 


( 8 ) 


u AB = 

0 


1 + TipVi + — r^rjpViVj H- 


e tot + Mi Vi + QgViV,- + 
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47re n 


■"tot 


v A (R) 

rfVi + QfAJiVj + 


m 


47ren 


e tot e tnt, + - ntufViVj + ^ itQ^YiYj + e^QjJ-ViV^ 


/<;'QjiY;Y,Y/.. - ii' i Q j \.Y,Y,Yi,. ) + Q'^Qj-iY : Y,Y h Y, + 


.Br^A 




/(*)• 


(9) 


Here e f, hf = S e 0 r i0i Qfj = efr^r^, and so on. Since atoms A and £> are 

£ P 0 

both neutral, = e^ t = 0. Accordingly, 

1 




47re n 


rtn fViVi-1 ) V.V,V t +Q5QgV.V,V t V,+ 


.AnB .BnA 


A^B-y 


/(«) (io) 


can be expressed as a sum of dipole-dipole (dd), dipole-quadrupole (dq), quadrupole- 
quadrupolc (qq), and other contributions. 

To evaluate the interaction energy between ground-state rare-gas atoms A and B , the 
classical interaction energy given by Eq. (jTOjl should be first converted into quantum mechan¬ 
ical operator. Perturbation theory may then be adopted to obtain the various perturbation 
contributions to the interaction energy at large R. 

Let the Hamiltonian of an isolated rare-gas atom X (.X = A, B) be H x . The Schrodinger 
equation 

H x tf = ( 11 ) 
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is solved for the n th excited-state energy E A and wavefunction ^, where the n — 0 case 
refers to the ground state. Accordingly, the full Hamiltonian of rare-gas atoms A and B can 
be expressed as 

H = H a + H b + U ab . (12) 

The interaction energy between ground-state rare-gas atoms A and B can be calculated as 

AE int = E 0 -(E A + E B ), (13) 


where Eq is the ground-state energy of H. 

To circumvent the need for solving the Schrodinger equation with Hamiltonian E [, E$ 
may be expressed in terms of {E A ^ A ]E B ^ B }^ based on perturbation theory 0, Q. 
Since atoms A and B are well-separated, an appropriate unperturbed Hamiltonian is the 
sum of the Hamiltonians of the isolated atoms A and H, 


H n = H J 


H 


B 


(14) 


Consequently, 


H = H 0 + U ABl 


where U AB given by Eq. (HU is the perturbation. 


(15) 


A. Zeroth-Order Theory 

At large R , the effects of electron exchange are insignificant. Ac¬ 
cordingly, for the n th excited state, an d E= E A + E B , where the isolated 

atoms A and B are described by quantum numbers r and s, respectively. For the ground 
state, Tq 0) = V’oVif and Eq' 1 = E A + E B . Correspondingly, A E int = E 0 — (E A + Eq) rs 
(E g 0 -*) — (E a + E B ) = (Eq + Eq) — (Eq + Eq) = 0. Therefore, to obtain a nonvanishing 
A Ei n t, it is necessary to go beyond the zeroth-order theory. 


B. First-Order Theory 

En ] = (T^IHabIT®). Since in the Coulomb case, the isolated rare-gas atom X (A" = A, 
B) has no permanent multipole moments in its nondegenerate ground state 0. presumably 
this holds true for the erf interaction with a sufficiently large c o or for the erfc interaction 
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with a sufficiently small u. Accordingly, the dipole terms are vanished = 

the quadrupole terms are vanished = and so on. Therefore, the first-order 

correction to the ground-state energy is 
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(16) 


Accordingly, AE int = E 0 - (E* + E*) « (^ 0) + £< 1} ) - (£ 0 A + £ 0 B ) = £< 1} = 0. Therefore, 
to obtain a nonvanishing AE int , it is also necessary to go beyond the first-order theory. 


C. Second-Order Theory 


En ■* = — X) £ Afll fo" ^ • The second-order correction to the ground-state energy is 
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which is always nonpositive. 

From Eq. (ITU]) , if only the dipole-dipole contribution is retained, we have 


r jdd 

U AB ~ 


47re. 




( 18 ) 


13 

















Accordingly, the second-order correction to the ground-state energy due to the dipole-dipole 
contribution is 


171(2 ),dd 
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In Eq. (fT9j) . the (r = 0, s ^ 0) and (r/0,s = 0) terms are excluded in the summation, due 
to the vanishing dipole terms, i.e., (i{jq \nf\'i{j^) = 0 (X = A, B ). 


For the erf interaction, f(R) = erf ^ R ^ . 
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( 20 ) 


Since e “ 


decays faster than polynomials when a is large, 
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_ erf( oj-R) ^ _ 3 RiRj) at large A. Accordingly, 
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1 


[erf(cni?)]" 


(<% — 3 RiRj)(5r 
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(47re 0 ) 2 R 6 

I ^ I Vtf > (Vf I /if I ^o) tyo \^\^r) Wo I $ I V>f) 


AE 

r^O s^O 


(21) 


B2> + E? - EH - E* 


• S), 


Similar to the Coulomb case (e.g., see Chapter 3 of Ref. [5j), we adopt the rotational 
average of {^|p^|^o l )(V , o 1 |A t ^|'0i 4 ) = \&a'\{' l P^\l iA W’o)\‘ 2 ^ and the rotational average of 
(V’f \vf iV’oWoVylV’f ) = KyKV’f |^W>| 2 , where = E e « r « and = E e JV 

a /? 

Also, note that 

3 3 3 3 3 3 

y y y y (e ~ 3 - 3 R i fRy)8 ii '5 jj > = y y (e - 3R i R j y= 6 . ( 22 ) 

Z=1 J=1 z' —1 j' — 1 1=1 J = 1 

Therefore, from Eq. (]2T]) . 

jpM,dd ^ 1 [erf(wR)] 2 ^ ^ |(y|p A |^)| 2 |(^f|p B |^)| 2 


247r 2 eg 


R 6 


EE 

r^O sj ^0 


B; 4 + b? - B 1 - B» 


(23) 


From Eq. (ITU]) , retaining also the dipole-quadrupole, quadrupole-quadrupole, and 
other contributions will produce additional terms in Eq. (ITT]) , involving V;VjVfc/(.R), 
VjVjVfcV;/(i?), and so on. For the erf interaction, it can be shown that V,;V ? /(A) 
decays more slowly than V*V jWkf(R), and ViV,-Vfc/(.R) decays more slowly than 
VjVjVfcV;/(-R), and so on. Accordingly, Eq 2) « E^’ dd at large A. Therefore, in the 
second-order theory, the interaction energy between rare-gas atoms A and R at large 
-R is 

= £ 0 - « + O ~ (4 0) + 4 1} + ) - (Eo + E*) = E ( 0 2) « 4 2) ’ dd 

1 [erf(o;R)] 2 |<y|p A |^)| 2 [<Vif|p B |^)l 2 (24) 

247 t 2 6 q # + ’ 

which has the [erf(cni?)] 2 R~ 6 asymptote. 


• For the erfc interaction, f(R) = crfc C' jR ) . 

In the second-order theory, the interaction energy between rare-gas atoms A and B , 
A E int = E 0 - (E* + E*) « 4 2) = - E K ^yyF ))|2 , is always nonpositive. 
Therefore, it is necessary to go beyond the second-order theory to describe the repulsive 
interaction energy at large R (as discussed in our paper), which is, however, beyond 
the scope of our discussion here. 
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FIG. 1. Interparticle interaction as a function of interparticle distance (in atomic units). 



FIG. 2. Potential energy curves of the He-He dimer associated with the long-range interparticle 
interactions erf (ur)/r, calculated using the corresponding CCSD(T). The ui = oo case is equivalent 
to the Coulomb interaction 1/r. 
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FIG. 3. Potential energy curves of the He-He dimer associated with the long-range interparticle 
interactions erf(cur)/r, calculated using the corresponding CCSD(T), CCSD, MP2, HF, PBE, and 


LDA. The uj = oo case is equivalent to the Coulomb interaction 1/r. 
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FIG. 4. Potential energy curves of the He-He dimer associated with the short-range interparticle 
interactions erfc (ur)/r, calculated using the corresponding CCSD(T). The ui = 0 case is equivalent 
to the Coulomb interaction 1/r. 
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FIG. 5. Potential energy curves of the He-He dimer associated with the short-range interparticle 
interactions erfc(o;r)/r, calculated using the corresponding CCSD(T), CCSD, MP2, HF, PBE, and 
LDA. The u = 0 case is equivalent to the Coulomb interaction 1/r. 



FIG. 6. Potential energy curves of the He-Ne dimer associated with the long-range interparticle 
interactions erf(cur)/r, calculated using the corresponding CCSD(T). The ui = oo case is equivalent 
to the Coulomb interaction 1/r. 
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FIG. 7. Potential energy curves of the He-Ne dimer associated with the short-range interparticle 
interactions erfc (ur)/r, calculated using the corresponding CCSD(T). The ui = 0 case is equivalent 
to the Coulomb interaction 1 jr. 



FIG. 8. Potential energy curves of the Ne-Ne dimer associated with the long-range interparticle 
interactions erf(wr)/r, calculated using the corresponding CCSD(T). The u> = oo case is equivalent 
to the Coulomb interaction 1/r. 
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FIG. 9. Potential energy curves of the Ne-Ne dimer associated with the short-range interparticle 
interactions erfc (ur)/r, calculated using the corresponding CCSD(T). The ui = 0 case is equivalent 
to the Coulomb interaction 1 jr. 



FIG. 10. Potential energy curves for the lowest triplet states of H 2 associated with the long-range 
interparticle interactions erf(wr)/r, calculated using the corresponding CCSD. The u = 00 case is 
equivalent to the Coulomb interaction 1/r. 
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FIG. 11. Potential energy curves for the lowest triplet states of H 2 associated with the short-range 
interparticle interactions erfc(o >r)/r, calculated using the corresponding CCSD. The uj = 0 case is 
equivalent to the Coulomb interaction 1/r. 















